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■ABSTRACT 

In this report our main concern is Image 
reconstruction of the object, emitting acoustic 
signals, using bispectral holography. Recording 
of the hologram using conventional techniques as well 
as using bispectral analysis is given. Conventional 
method is employed for image reconstruction from the 
recovered hologram signal. The bispectral holographic 
imaging technique has the potential theoretical 
advantages of ideal recovery in the presence of noise 
and distorting medium. Computer simulation of this 
processing strategy carried out in this work confirms 


this. 



CHAPTER 1 


INTRODUCTION 

Holography is a technique, meant for the purpose of 
recording and displaying an image in its true three-dimensional 
format. It is a synthesis of two branches of optics, namely, 
interferometry used in recording the hologram, and diffraction 
used to display the image. This technique of photography has 
come to be known as holography after the word 'hologram' coined 
by Gabor to describe the inter ferogram that is used to record 
the wave front to be reconstructed. The initiation for this 
synthesis came from engineers and physicists working in the 
area of communications. In communications the analogies to 
interference and diffraction are respectively modulation and 
demodulation. It is obvious from this analogy that the 
information stored in an interference pattern (modulation) can 
be retreived as a result of diffraction of light by the 
recorded interference pattern (demodulation). 

In the following section a simplified discussion of 
holography based on interference and diffraction is undertaken 
At first, the interference phenomena is studied. As shorn in 
Pig. 1.1, a plane monochromatic wave is split into two equal 
amplitude plane waves which then traverse separate paths 
until they are recombined. If the mirrors and beam splitters 
are perfect and if the optical path lengths are equal, the 
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output beam will be exactly like the input beam except for 
the loss in beam encountered at the beam splitter. A screen 
placed in the beam will therefore show uniform intensity. 

If one of the mirrors is not quite parallel to the 
wavefront due to a misalignment, the reflected wave from 
it will combine with the other beam to form a linear system 
of fringes as shown in Pig. 1.2. Destructive interference 
occurs who rover the path length difference between the 
wavefronts is an odd multiple of half the wavelength. 

If one of the mirrors is replaced with a concave 
parabolic reflector, the interference pattern wo would 
expect to see is shown in Pig. l.J. Fote that this ring 
pattern shows decreasing spacing with increasing distance 
from the centre. This type of pattern is called Fresnel 
Zone pattern. If we impart a tilt to the spherical wave, 
we obtain a partial Fresnel Zone pattern as shown in 
Fig. 1.4. 

When an aperture of some kind is interposed in a 
collimated beam of light, one might expect to see on the 
screen, a pattern of light sharply defined by the shadow 
of the aperture. Looking closely, however, the observer 
notes that some light exists in the shadow zone. The study 
of this phenomena is known as diffraction theory. The 
os sense of the theory known as Huygens 1 principle is that 
each point on a wavefront can be considered as an elementary 
point radiating a spherical wavefront. When the elementary 
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wavelets are summoned' in the prescribed manner the result 
is identical to the actual wavefront farther down stream. 

This principle is useful in analyzing the behaviour of the 
light in the presence of various types of apertures and 
obstructions. 

A film containing the binary interference pattern 
generated in Fig. 1.2 is placed in the path of a collimated 
monochromatic light beam as shown in Fig. 1 . 5 * In this 
figure we consider a single wavelet in the centre of each 
transparent space. Note that each circle represents an 
equiphase surface advanced in phase by 2a from its neighbour 
to the left. The spatial separation of each surface from 
the next is one wavelength X . We can draw a plane phase 
front tangent to the elementary wavelet phasefronts , three 
of which are shown in Pig. 1 . 5 * For the binary pattern x^e 
can draw many more plane wavefronts at ever increasing angles. 
Comparing Fig. 1.5 with Pig. 1 . 2 ? we see that we have 
reconstructed the beams that made the interference pattern 
in tho first place, with one exception? we have an extra 
beam at the conjugate angle. 

The second example of the interference pattern was 
the Fresnel Zone pattern (Fig. 1 . 3 )* Fig. 1.6 shows what 
happens when this is illuminated. The Fresnel Zone pattern 
is, in effect, a locally linear grating, acting much the 
same way as our first example, except that the spacing varies 
as described previously. Hence the locally diffracted light 
changes in direction and the total effect provides a spherical 
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■5 RECONSTRUCTION OF THE PLANE WAVES' 
THAT MADE THE INTERFERENCE PATTERN 
IN FIG. 1 -2 
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wavefront, again duplicating the beams that caused the 
original interference pattern. Again we have an extra 
wave of opposite curvature. 

The last example considered was that of an offset 
Fresnel Zone pattern (Fig. 1.4). For completeness we show 
in Fig. 1.7? the result of introducing this pattern in a 
beam of light, although one can anticipate the consequences# 

The significance of the foregoing discussion is 
that interferometry has always had the capacity of recording 
a wavefront and later reconstructing it. These three 
examples chosen above clearly explain how holography works# 

The first example (Fig. 1.2) was intended to demonstrate 
both interference and diffraction. The second example 
(Fig. 1.3) > a slightly more complicated situation, was 
chosen to show how Gabor Holography works and why it is 
not successful. Using the Gabor approach, one is recons- 
tructing three beams of light, all occupying the same volume. 
That is, if one considers Fig. l.S, one is forced to look 
at the plane wave and the converging extraneous wave both 
of which tend to mask the light from the desired wave, when 
looking for the replica of the spherical wave used in making 
the interferogram. 

The third example, Fig. 1#4» shows how this difficulty 
is overcome hy using an offset in the spherical beam. The 
same result can be produced by bringing the plane reference 
wave in at an angle, although in this case the interference 
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pattern is not a circular Fresnel pattern. In this case, all 
the three diffracted beams are separated in space , thus 
completely solving the problem of overlap. 

To make these discussions conform to what one generally 
likes to think of as holography, one further step is required? 
namely, to record and display an image, of a complicated object. 
This step only requires that we accept Huygens' principle which 
states that any arbitrary wavefront representing the light 
reflected from a complicated object can be considered to be the 
sum of a large number of point sources distributed over that ■ 
wavefront. Since we have shown that we can imago a single 
point or a spherical wave, it is reasonable to suppose that 
a distribution of points can also be imaged. 

The foregoing theory depends upon the overriding 
principle; that is, that the two beams must be capable of 
interfering. This, in turn, requires a property of light, 

9 

known as coherence. A wave is said to be perfectly coherent 
if its phase structure is fixed with respect to time. The 
degree of coherence of wave influences the recording of 
interference fringes and consequently the formation of holo- 
gram. In recording a wavefront, the sensitivity or time 
constant of the detector determines the integration time 
necessary for recording of interference fringes. If the 
relative phase of the wave changes during the recording time, 
the fringes move and quality of the recording is degraded. 

If the phase of the wavefront fluctuates randomly, the 
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recording of fringes is best described in terms of the mutual 
coherence function which is defined as the cross-corre lation 
of the complex functions describing the wave at two points 
in space at different times under the assumption that waves 
are all polarised in the same sense and stationary in time. 
Lasers are commonly thought of as single frequency, highly 
coherent optical sources. However, they are often multiple 
frequency sources with relatively low coherence. By proper 
precautions and controls, a laser with a temporal and spatial 
coherence properties can be obtained. Sources for Microwave 
and acoustical wavelengths are normally coherent. 

bo far we have discussed optical holography but our 
main interest in this thesis is in some aspects of acoustical 
holography. Hirst, we will see fabat are acoustic signals. 

The term 'Acoustic' applies to all mechanical vibrations and 
compressional and shear waves having any frequency. A very 
broad spectral range is encompassed within the field of 
acoustics. Any mechanical vibration with frequencies ranging 

1 p 

from a lower limit of almost 0 Hz to an upper limit of 10 c Hz 
is within the realm of consideration. 

How we will relate the acoustical holography to optical 
holography. It should be noted that the interference and 
dif faction are phenomena common to all forms of energy whose 
propagation can be described in terms of a wave motion. All 
that is required is the ability to record an interference 
pattern in the particular type of radiation we happen to be 
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using and. use this pattern to diffract light into a replica 
of the original wave. In this way it is possible to produce 
a visual image of an object irradiated with invisible energy. 

Photographic film is insensitive to acoustic radiation, 
Point by point sampling method using one or more acoustical 
receivers a.nd liquid surface holography are some of the methods 
developed for recording the acoustic hologram. The recording 
procedure of these methods is explained in detail in the 
second chapter. 

Thus we have studied the recording of the hologram 
and reconstruction of the original wave (object wave) on a 
qualitative physical basis in this Chapter. Chapter 2 deals 
with the hologram recording and image reconstruction, in 
mathematical terms. The third chapter discusses the basic 
principle behind the bispectral passive imaging and the 
computational aspects involved in estimating the bispectral 
density. The fourth chapter presents the simulation work 
that lias been done for the reconstruction of images of some 
idealized objects. The fifth and final chapter discusses 
the computational results and recommendations for future work. 
The computer programmes are listed in the Appendix. 
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CHAPTER 2 

C OK VE NT I OH AL HOLOGRAPHY 

In Chapter 1, how a modified form of interferometry 
could he used to record a wavefront wa. s discussed on a quali- 
tative physical basis* This record, in turn could he used 
to reconstruct the wavefront by the process of diffraction. 
This method of photography has come to he known as holography 
after the word 'hologram' coined by Gabor to describe such an 
interferogram. In this chapter the process is described in 
mathematical terms. 

2 . 1 Recording Process for Optical Bigna.ls 

Consider two beams of radiation intersecting on a 
detecting plane as shown in Pig. 2.1. These two waves come 
from two transducers driven by the same oscillator. We can 
describe the optical field at a point (x,y) on the detecting 
surface due to beam 1 as 

Sl(x,y) = a. 1 (x,y) cos|>t + 9 1 (x,y)} (2.1) 

and the field at (x,y) due to beam 2 as 

& 2 (x s y) = a 2 (x,y)cos[a)t + 9 2 (x,y)3 (2.2) 

where a(x,y) is the amplitude of the wave and <p(x,y) is 
the phase of the wave. The total amplitude at (x,y) is the 
sum of these two individual beams. The result is 

S(x*y ) - S 1 (x,y) + & 2 (x„y) 

= a 1 (x,y)cosCw t+ 9 - L (x,y)]+a 2 (x > y )cos[ot+ 9 2 (x,y )} 

( 2 . 3 ,) 
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Since no detector exists which can detect the oscillations 
at optical frequencies, the test that can be done is to measure 
the intensity. The signal that is actually recorded is therefore 
some function of intensity, which can be written 

I(x,y) = < [S(x,y)3 2 >£ (2.4) 

where < denotes a time average over many cycles of 

radiation. 

Substituting eq.(2.3) into eq.(2.4) yields 
i( x >y) — + &2 t "*■ (Sq&2^_ 

w he re 

i aj(x,y) < [1 + cos 2(ot + y ( x * y ) ) 3 ^ 

Sg = i a|(x,y) < [1 + cos 2 ( w t + <p 2 (x,y))] > t 

(2.5) 

(^l^f* i a 1 (x,y) a 2 (x,y) < cos(2wt + <p 1 (x,y) 

+ 9 2 ( Xj y)) > t 

(SiS 2 )„= i a x (x,y) a 2 (x,y) <co s(<jJ 1 (x,y) - <p 2 (x,y)) > t 

This equation clearly shows why the temporal frequency 
term is so easily dropped. 1 An oscillating function such as 
cos(w t) averages to zero, if many cycles are considered. 

He nee , e q . ( 2 . 5 ) re duce s to 

I (x,y ) = -H[a£(x,y) + a 2 (x,y) + a- L (x,y) a 2 (x,y) cos^Cx^) 

- <p 2 ( x »y)) 3 } (2*6) 

Note that we have succeeded in preserving phase 
information even though the recording was done on a pha.se 
insensitive medium. This, of course, is the secret of 
holography. 
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2.2 Reconstructio n Pro cess for Optical Sig nals 

In Chapter 1 , it is noted that the wavefront stored 
on the interferogram could he reconstructed hy the process 
of diffraction of a plane wave. This plane wave' can he 
expressed as 

= a^(x,y) cos( at + c p^(x,y)) (2.7) 

Suppose we have exposed a photographic plate to the 
intensity shown in eq.(2.6). After development it possesses 
a transmittance proportional to the intensity. Then the 
wave after passing through the plate is modified hy it as 

KI(x,y) 8 j(x,y) = -[ a£(x,y) + a|(x,y)3 a^(x,y )cos[o t-h(p^(x,y )] 

+ - a 1 (x,y) a 2 (x,y) Ccos[wt+cp 1 (x J y)-(p 2 (x,y) 

+ 93(x,y)] + cos[ ut-<p 1 (x,y)+9 2 (x,y )+(pj(x,y )] } 

( 2 . 8 ) 

where IC is a, constant. 

If we consider the terms one at a time we note that 
the first term represents the reconstruction wave modified in 
amplitude. The second term assuming that 8^ is a replica of 
8 2 f he c ome s 

4 a|(x,y) a 1 (x»y) cos[ w t + 9 1 (x,y)3 

which we recognise as the original wave, 8^, modified in 
amplitude. If S p h^pgns to he a plane wave, a p is a constant 
and a perfect reconstruction is obtained. With the same 
reconstruction wave the third term of eq, ( 2.8) he comes 



19 


f a 2 ( x >y) a i( x jy) Q°s[ wt - <p 1 (x J y) + 2g> 2 (x,y)] 

This is the extraneous term we referred to earlier. 

There are several more options available. For example, we 
could let he a replica of S 2 with opposite curvature,’ 
that is 

S*(x»y) = a 2 (x,y) cos[ ut - <p 2 (x,y)] 

Then the second and third terms of eq.(2-8) become 

ag(x,y) a 1 (x,y) cos[ ot + cp 1 (x,y) - 2rp 2 (x,y)3 
and 

- a|(x,y) a 1 (x,y) cos[ wt - ^(xjy)] 

Here we see that the last term becomes a re construction 
of with opposite curvature and the second term is the 
extraneous imago. Another possible form of Bj is giving 
reconstructed terms. 

•| a£(x,y) a 2 (x,y) co s [ ut + 2<p 1 (x,y) - <p 2 (x,y)] 

and 

| a l( x j y ) a 2 (x,y) cos[o t + <p 2 (x,y)] 

Thus wo have reconstructed S 2 and demonstrated the 
complete ambivalence of the interferometer. While recording 
the wavefront of one arm wo have also recorded the other. 

Hence either : one may he re constructed by diffra cting the othe r 
beam through the interforogram. Additionally by using the 
conjugate of the one of the interfering beams as the diffracting 
beam we can re const ruet; the conjugate of the other. 
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2 * 3 Re cording of an Acoustic Hologram 

In case of acoustic signals, because of their long 
wavelength, measuring the phase is not a difficult problem. 
Detectors capable of measuring the amplitude of oscillations 
of the acoustic signals are readily available. The system 
used, for recording the hologram is shown in Fig. 2.2. Here 
the received signal (scanned signal) is multiplied by the 
oscillator signal (reference signal) and sent through a low- 
pass filter to yield the desired hologram signal. The 
mathematical procedure is as follows: 


S- L (x,y) - a 1 (x,y) cos[ ut + cp 1 (x,y)l (2.9) 

After multiplication by a 2 cos(wt + c p 2 )> this becomes 
a ? a 1 (x,y) c 

8’(x,y) 2 i c os [ 2 wt + cp 1 (x,y) + <p 2 J 

+ cos[cp 1 (x,y) - 9 2 3 J (2.10) 

After low pass filtering, we are left with 


&" (x, y ) 


a 2 a p(x,y) 

2 


cos[ 9 1 (x,y) - 9 2 3 


( 2 . 11 ) 


Comparing eqs.(2.9) and (2.11) we see that we have 
succeeded in obtaining the same information without the 
attendant extraneous terms a| and a^. 

The signal ;S" (x,y) may be stored on a magnetic tape. 
The only way to get the full amount of information from this 
data is to reconstruct the wavefront. To do this, the data 
must be recorded in the form of a hologram. The basic method 
is to intensity modulate a light source with the stored signal 
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on magnetic tape about some average intensity, move it in 
position and velocity in proportion to that of the detector 
and record the light on a light-sensitive medium that is 
shaped proportionately to the surface over which the receiver 
was scanned. 

In practice only planar surfaces have been used, but 
scan lines in the surfaces have been straight lines or circles. 
The most common light sources have been an incandescent bulb 
focussed to a very small point, and the face of a cathode ray 
tube with a moving intensity modulated spot. The first method 
requires a mechanical scanning system whereas the second one 
requires an electronic scanning system. Each method requires 
a camera to photograph the moving light source. Reconstruction 
process is same as the one explained for optical signals. 

2.4 ' Liquid surface holography 1 

'Liquid surface holography' is another means of recording 
the interference pattern. The system used for imaging is shown 
in Fig. 2.3. It includes an object beam transducer that 
insonifies the object uniformly. The object scatters, diffracts, 
absorbs energy, and modifies the phase distribution in the wave- 
front. A pair of acoustic lenses are used to image the sound 
distribution of the object into the plane of the hologram. The 
geometry of the wavefront of the reference beam should be simple 
so that the wavefront oan be easily simulated in a beam of 
light. This means that the wavefront should be either plane, 
or spherical. The object beam and the reference beam should be 
essentially coherent. The rapid response of the liquid surface 
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permits a few percent difference in frequency between the 
object and the reference beams. The sound field of the 
object beam mixes with the sound field of the reference beam 
transducer to form an interference pattern at liquid-air 
interfa.ee in the imaging tank. The interference pattern is 
imprinted in the liquid surface as variations in elevation. 

An imaging tank is used to isolate the hologram from 
liquid surface disturbances in the main (object) tank. Use 
of a imaging tank also allows the use of different liquids 
for hologram formation and for object immersion. Thus? 
liquids can be chosen which are optimum for each of these 
purp oses. 

The hologram surface imposes phase variations in the 
reflected beam of light which cause diffraction of the light. 
First order diffracted light is passed at the spatial filter 
which blocks all other light reflected from the hologram 
surface . 

Lens 1 is used to image the point source of light at 
the plane of the spatial filter after reflection from 
hologram surface. Lons 2 images the hologram surface on the 
ground glass viewing screen. Since the hologram is a, focussed 
image type, the ground glass screen displays the desired image 
of the object. Although the hologram modifies only the phase 
of the light, there is amplitude imaging because the amount of 
light diffracted into the first order varies as the amplitude 
of the interference pattern on the liquid surface? which 
inturn, varies as the amplitude of the wave in the object beam. 
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2.5 Recording of the Hologram Employing Digital Signal 
Processing Methods 

An alternative solution for the reconstruction of 
the wavefront of the object wave is using digital signal 
processing methods. Here ? we use two receivers; one is 
fixed and the other scans over a hologram plane. 

The signal detected by the scanning receiver can be 
expressed as 

S 1 (t,x,y) = A(x,y) cos[ at + cpx(x,y)3 + n-^t) 

where (x,y) is a point on the hologram plane and nj(t) is 
the additive Gaussian noise. A (x,y) and <p(x,y) are amplitude 
and phase distributions of the scanned signal. 

The signal detected by the fixed receiver 

S 2 (t) = reference signal = B cos( wt + c p 2 ) + n 2 (t) 

where n 2 (t) is additive Gaussian noise. Both n^(t) and n 2 (t) 
are assumed to be independent of the received signals * but not 
so in between them. .1 ?■' - 

The cross-power spectral density between the object 
and reference waves is computed at all the scanned points on 
the hologram plane and is given by the following expression; 

Cg g ( to, x»y ) = A(x,y) B cos(tp 1 {x,y) - <P 2 ) + c n n 

where C -is cross-power spectral density between the 
n i n 2 

noises ( t ) and n 2 (t). 

Xt is clear from the above expression that the cross- 
power spectral density C g Q contains the amplitude and phase 
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information of the object wave. Hence the cross-power 
spectral density computed at each point on the hologram 
plane forms the hologram signal in the Fourier domain. The 
image is retrieved as a computer print- out, simply by 
performing an inverse Fourier transform on the hologram 
signal. 

The signal processing carried out in the computer 
is shown in Fig. 2. 4« 
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^ C- 


<L^ 


B I SPEC TRAP HOLOGRAPHY 


The amplitude and phase of the object wave can 
be obtained by cross correlating the object wave and 
reference wave on the hologram plane. This is the 
conventional way of recording the hologram, which was 
discussed in the last chapter. If the signals illuminating 
the object or object emitted signals are non-Gaussian the 
phase of the signals can be derived from bispectra or 
higher order spectral analysis. The condition that signal 
should be ' non- Gaussian is satisfied in many practical 
applications, since any periodic signal can be regarded 
as a non- Gaussian signal and self-emitting signals from 
objects such as complicated mechanical systems can also be 
considered as non- Gaussian signals. 

Two_re ceivers one fixed and the other scanning 
over the hologram plane are used to obtain the holographic 
information. Auto-bispe ctral and cross-bispe ctral analysis 
are carried out on the detected signals. By taking their 
ratio the required information is obtained. This is 
explained in the later parts of this chapter. Image 
reconstruction is performed by using conventional holographic 
methods. 

Bispe ctral analysis provides a means of obtaining 
holograms free from additive Gaussian noises in the amplitude 
and shape of the power spectra. 
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The application of bispectral analysis is equivalent 
to the utilization of mutual relations among three frequency 
components, that satisfy f^+fg+f^ = 0. In the case of a 
Gaussian signal each frequency component is independent. 

These facts produce very effective and special features 
when this method is applied to holographic imaging systems. 

The principle of bispectral analysis for getting the 
hologram signal is dealt with in the following sequence. 

First the descrpition of the data acquisition system and 
assumptions postulated for generalized object illuminating 
signals and additive noise are given. Secondly by discussing 
the properties of the bispectra of the signals detected at the 
hologram plane, the procedure required to obtain the hologram 
signal is derived. Finally the image reconstruction process 
is discussed briefly. 

3.1 Theory of Bispectral Imagingf 23 

Consider the sj^stem shown in Fig. 3*1^ where non- 
GaussianrMid om signals S i (p,t) (i = are assumed to 

be emitted from the object plane n^. n(t) is the noise added 
in the channel . These signals are received at the hologram 
plane by two receivers. One. is a fixed receiver R-p> 
which is assumed to be located at the origin and the other 
one is scanning receiver R , located at r 2 on the same plane. 
n 1 ( t ) and n 2 (t ) are the noises added to the dete cted signals 
at receivers lip and R 2 ^respectively. The signals on the 








plane may "be considered as sum of the following three 
comp onents . 

a) Signals from the object plane. 

b) Surrounding noises n (t), n(t). 

s 

c) The additive noises n^(t) s ( "t ) which are added 
at the corresponding receivers. 

Information about the amplitude and phase of the 
wavefront of a given frequency at each point on the hologram 
plane is required to construct the hologram. 


Assumptions ; The assumptions postulated for the 

object illuminating signals and the transmission medium 
arc given be lows 


i) The object illuminating signals Sj_(p*t) (i=l,..»M) 
are harmonizable (this means that ( p>t) can be expressed 
as the first expression of Eq.(3.1)) and third order 
weakly stationary non~Gaussian signals and without loss of 
generality they are assumed to have zero mean. 


S.(p,t) = / exp( j u) t) dZ ± ( w) 
■*” 00 

< azpw) > = o 


(3-1) 


W *2> d2 “ = <dz h“l) aZ i ( “‘3 ) j +B +B =0 

12 3 

who re < . > d en ote s an ensemble averaging ope rati on and 
b . ( o w ) represents the bi spectral density of) S i (p,t);..? 
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ii ) Noises are assumed to be Gaussian. Eo restriction 
is made about their power spectra. It is also assumed that 
they are statistically independent of the object illuminating 
signals Sj_(t) (i = l 9 ...M). Eo assumption is made about the 
independence among these noises. 

The received signals u (r , t) (q = 1,2) at the two 

tL 4. 

receivers can be expressed as follows : 

00 M . . 

^(r q9 t) = / / l A i (p)Z i (w )H(w } p,r q ) e 3u) ' c d m dp 

TC-, - “ 

+ n(t) + n (r^jt)' , q=l,2 (3*2) 

where A(p) is the amplitude of the object wave at a point p 
on the object plane and H( cd ? p, r) is given by 


oo 

/ h(x ? p, r)e” 3wi: dx 

co 

where x is the time instant at which you apply the input and 
r is an observation point on the hologram plane. 

Under the assumption that turbulences of the medium 
affect only as random pha.se shifters on propagating waves 
and their variation in time are so slow that they can be 
regarded as fixed values while -bi spectral' analysis for each 
data is carried out, the transfer function H( to ; p,r) is 
expressed as follows; 

G( CD ) r , . . 

H( w; p, r) exp[-jkd(p,r)jexp[ ja(u; p,r)J 

d X d(p,r) , 
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where x , k denote respectively, the wavelength and wave- 
number of a wave with angular frequency o> ; a( w j p,r) is 
the accumulated phase shift due to random medium between 
p and r ; G(w) is the amplitude transmitting factor. 

Another assumption made on channel is the accumulated 
phase shifts caused by random medium are subjected to the 
same probability distribution for any p and r, so that they 
are expressed, as a random variable with index to as b(«). 


The transfer function of the medium described above 
does not take the multipath effect into consideration. The 
reason for this is that the size of the random sc^tterers 
is much greater than a wavelength. Because of that the 
wave scattered by the particle is largely confined to a 
small angle in the forward direction. 


V r < 


<T 


Substituting H(u ; p,r) inBq.(3. 
M 00 

t] = iii l [ 4 na(P.i- q ) A i (p) 


2), we get 

exp[-jkd(p,r q )3 


exp[ja{(o;p,r q )3dp] exp(j«t)dZ i ( co)+N q (r q ,t) 
q = 1,2 (3-3) 

where N (r ,t) (q = 1,2) are the overall corrupting noises 
at the receivers and they are independent of S^(p,t) 

(i = 1, ... M) from the assumption (ii). A(p) is the 
amplitude of the emitting signal at a point p on the object 
plane * 

Consider the signals X^(t) = X(t) + Yj_(t) (i = 1?2) 
whe re Y. ( t ) ( i — 1 , 2) are weakly.* stationary Gaussian signals 
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and X(t) is a third order stationary non-Gaussian signal 
and independent of Y^(t) (i = 1,2). Then the third order 
cross-correlation function R v v v ( x-, s r ) between X-, ( t ) 
and X 2 (t) is given as follows; 

Rl l X 2 X l^ Ti ’ T 2' > “ < ^l^)" <2 1 ( t ) > 3[2 2 (b+T 1 ) - < X 2 (t+'c 1 ) > ] 

[X- L (t+'u 2 ) - < X^( t+Tp) >]> 


" r xxx^ t i ,t 2^ + a? 1 Y 2 Y :J T i ,T2 ^ 


R XXX^ t 1 ,x 2^ 

where the fact that the third order correlation function 
of any Gaussian signal is identically zero is used. 

By applying this relation to u (r^t) (q = 1,2) in 
4>q.(3*3), the third order cross- correlation function 
R ( t, , t ) which is abbreviated as R-, oi ( T i 9 x o ) can de 

\3.~| _L 2 Jm J- Cl 

derived as 


M oo 


ipfa-v) = h J_, .. L L a “i l a “2 


1 2l v 1 f 2 7 1^=! i 2 k i 3 =l ioo 


Shii j 
2 0 


sc.p ac 2 ) 


5 'IT -1,-1 ±1 2 2 

~y x 1 x 2 (x x + \ 2 ) 

<^exp 3 [b(w^)+h( tog ) +b ( *— to-j_ 

z Ai^pp ■ A d (p 2 : i 

*1 — T exp(-jk 1 d 1 )ap 1 n. — ‘ e x p(- 0 k 2 a 2 )a 


*1 a. 


whe re 


/ — 2L exp(j (k 2 _+k 2 )d^)dp^ 

1 ■ ' 3 ; TfvT 

d (Pl» r 2^ and d i " ■ d.-CPi» r i ) ( i= 2# 3) 


(3.4) 
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Then the cross-bispe ctral density ’ whic:i1 is defined 

as the Fourier transform of &]_2l( T l' ? can ex P resse( ^ as 


M K 


a(« 1 )G(o 2 )G(-a) 1 -aj ? ) 


( V ») 2 I I r v ± , « (. lf ._) 

i-i = l ±w=l i„=l 1 2 -i A ^( , XT- L + X ^ 1 ) 1 


‘1 2~ 3" 


'1 2 1 


■J-V 2 ( X 1 + *2 ) 


< exp|j[b((o 1 ) + b(o3 2 ) + b(-a) 1 - “ 2 )3\> / — ii£l ) exp(-jk 1 d^)dp 1 


TTj u l 


r A i 2 ( P p ) c •^■ix(P^) 

J — exp(-jk 2 d 2 )dp 2 / —g — exp(j (k x +k 2 )d 5 )dp 3 

12 TC t *3 


”l 3 


(3-5) 


In the same way, the auto-bispe ctral density W 2^ 

obtained from hq.(3*5) and is given belows 


B lll(“l’" 2 ) 

0 M M M G(O) )G(to )G(_a) _ 0 ) )' 

= (l/2n) Z l . I J b i i (<V“ 2 > Y 

1 1“ 1 1 2~ 1 i 3 _1 12 1 -j X i X 2 X 1 +X 2 ) ;; 

< exp {j[b(« 1 )+b(« 2 )+b(-» 1 -« 2 )^> / e xp (-jk 1 a 1 )dp 1 


AjgCPg) 


e X p(-jk 2 d 2 )dp 2 / xp ( j ( k x +k g ) d dp ^ 


*1 2 


(3.6) 


Equations (3*5) and (3.6) give the bi spectra of the signals 
at the hologram plane and they are expressed as the product 
of the following three main terms ; (i) the term which 
depends on the bispectrum of the object illuminating signals 
S i (p f t) (i = 1,.. . M), (ii) the term due to the phase 
turbulence at three angular frequency components {0 1 ^ a) 2 sC °3 
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■where + u ^ = 0 ; rind (iii) the object dependent term. 

In the expressions (3-5) and (3*6) the difference is only the 
use of or d^i thus if the r n tio between them is taken? the 
o the r fa c t or s ge t ca nee lie d . 


D erivation of Hologram Signal from Bi spectral Analysis : The 
data required for hologram signals are the relative amplitude 
and phase of the object waves at the hologram plane. They may 
be expressed as follows: 

A g (w;r) exp[j<p s (u>;r 
A r (w;r) exp[ jcp r (w ; r )3 


A 0 ( w ;r ) exp[ j<p Q (w ;r)3 


A ( oj ; r) 

exp 

A r (“ } r) 


j[9 s ( a) ; r ) “ 9 r ( w ? r )3 

(3-7) 


where A r exp(jcp r ) is a reference wave and A g exp(jcp s ) is the 
object wave at an observing or sennning point r on the plane. 
Eqs.(3.5) and (3.6) show that the hologram signal at a 
scanning point r^ can be immediately given by the following 


ratio: 


WvV 

B u:lll(“l’ “2> 


M A ii (P t ) 

l / ' 1 — exp(-jk 1 d 1 )dp- 

i l =1 *1 1 


M Ap-, (p-, ) 

I / 5~™ e xp(-j k 1 d 1 ) dp-L 


1 1=1 *1 


A o( w l> exp C i 9o( t * ) i» r 2 ^ (3-8) 


where B u ,i L -^ 1 ( u ) is assumed to be non-zero, 
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Equation (3«8) implies that the hologram signal can 


be obtained as the ratio of the or os s-bi spectrum and auto- 
bispectrum of signals at two points on the hologram plane. 

In this expression, only the desired factors are obtained 
and the terms which depend on the bispectrum of the signal 
and the turbulences are cancelled because the ratio is used. 
Equation (3*8) also indicates that the result is free from 


additive Gaussian noises. 


Image is reconstructed by carrying out the following 
operation: 


I(P') 


A (w,; r) 

/ e xp[ 39 0 ( w i»r)3 exp[jk d(r,p ' )3dr 

jx 2 d(r*p’) 


where p* denotes an image point to be reconstructs! and 
d(r,p') denotes the distance between r and p’. 


3*2 The Computational Procedure Involved in Estimating 
t he Bispectral Density 


The signals detected by the scanning receiver and 
fixed receivers are x^(t) (object signal) and Xp(t) (reference 
signal )re spe ctively. Their third order correlation functions 
of Xj(t) and Xp(t) are assumed to be ergodic. 


Then the third order auto- correlation function x^(t) 
(q - 1,2) is defined as 

1 T 

"^lll ( T l' T 2 ) = Lim T / x q( t ) x q( t+^i )x q ( t+Tp ) dt 
T o 


q. = 1,2 


Hence observing the signal x^(t) (q = 1,2) over a period T, 
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one c^n compute 


R 111T (t 1 ,X 2 ) ~ T 


T 

/ 

o H 


X q (t + x ± ) x q (t + t 2 ) dt 

q = 1,2 


To compute %]_p( ' one Ras °^ serve cig na l 

over many such, intervals of time period T? and has to 


compute E -mp( T i > T 2 ) each time. Then the average of such 
r 111T< i: 1’' c 2> ,b gives an estimate of Rppp( , % 2 ) » whose 
Fourier transform is nothing but an estimate of the auto- 
bispectral density of x(t). An alternate way of doing this 
is to compute Bpppp(o}p , oi ) (Fourier transform of Rp^lT^ " r l i ' T p) ) 
each time after observing the signal over a period T and then 
taking averages of such Bpppp( ^p, w 2 ) ,s . Birst we will see 
what is Bppp T ( w p» w 2 ) in "terms of X( to ) . 

B U]t ( «p, w 2 ) is defined as R(R 111T ( * 1 * T 2 )) where 
F( . ) denotes the Fourier transform. Therefore, 


1 t 00 00 — jui- 1 T- 1 “3 10 0^0 

B 111T^ “l' w 2 ) = T f f f x(t)x(t+Tp)x(t+'t 2 )e B e ’ 


0 —00 —00 


T 


ji f f J x(±)x(x- L )e^ a>t x(r 2 )e 


dt dTp dT 2 


-j“ 2 ^ -inn 

0 


— 00 — 00 


e ~3 “2 T 2 dt dtp dTg 


m / x ( T l) e 

T -00 x 


dv. 


“ “3 VV 

1 dTp / x(T 2 )e 2 2 

mm OO 

T 3 (op+a) p )t 

/ x(t) e ^ ^ dt 

0 

[since x(t) is observed only for an interval T, 
outside this interval its value is zero. Hence 
we can take the limits from -00 to +od ] 


f X ( 0)p ) x ( » 2 )x*(w 1 +u 2 ) 
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The cross-bispectral density of x-j_(t) end x 2 (t) is 
defined as 

To get an estimate of ai^) , the some procedure 

that is discussed above for estimating w ) is ada P' fced 

here. Before actually going to the discussion of the computa- 
tional procedure involved we will express ®i2iT^ u l > ^ e:rins 

of xq(w) and x 2 ( w )* 

is definGd as ■ B '^l2lT^ T l' i: 2- ) - 1 * 


T » 


B 121 T^ W 1 ,W 2 ^ = f / / /x 1 (t)x 2 (t+T 1 )x 1 ( t+r 2 )e 


-3 w i T i “ , 3 a) 2' c 2 


dt dr-q dx 2 


1 T 00 

f / Xi(t) / Xp( t+r 1 )e 

o — CO 


'^ a) l T l 


“’*“2 T 2 


d^i / x 1 (t+T 2 )e * dr 2 


= | / ^(t) / m X 2 (t 1 ) * e 1 1 dT] 


I *i (x ? ) 


j 0) o t -;jiD 0 T 


2 2 


dr. 


T 


/ x i ( t )e 


j(«i+» 2 )t . s ;3Vi 


/ x 2 ( T l) e dT ] 


00 

"“3 

J x-i (T„)e * dr. 


fh - 2 } “ 2 

~ X 2 (oji) x 1 (u) 2 ) x 1 *(m 1 +a) 2 ; 


(3*23) 


Thus the expressions for 2 ) ’ B l2lT^ u l ,w 2^ are der:i - ved * 

which are the triple products of the Fourier coefficients of 
x . (t) (q = 1,2)* Hence to get good estimates of and 

B l2dPl ,w 2) 0n3 kas to average over many such triple products. 
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The computations may be surprisingly expensive despite the PFT 
Moreover? all the pitfalls .known from ordinary power spectrum 
analysis occur here too, some of them with new twists. Thus 
it seems worthwhile to describe and evaluate the different 
procedures in detail. 

Assume that the bi spectrum is to be estimated for 
frequencies between 0 and w 0 /2 with a bandwidth A 0 * Hence 
a sampling rate of w samples per unit time is atleast needed 
and uninterrupted records whose size is atleast « e Q /& Q > 
Prefiltoring of the frequencies above id 0 /2 is necessary to 
avoid aliasing. A total of K records (not necessarily 
consecutive), each of length U = MET 0 is assumed to be 
available. The value of N should be convenient for the FFT, 
say a power of 2 and M preferably is an odd integer (M — 2L+1) 
To satisfy both conditions a compromising value of N Q may 
have to be chosen. Also, it may be necessary to adjust the 
records, either by adding zeroes at the ends or by allowing 
some overlap between consecutive records, to reach the exact 
length N; the final results will have to be corrected for 
these adjustments. The total number of points is of the 

order ET = KN. 

tot 

Because of its symmetry, it suffices to estimate 
the: bi spectral density B(w 1 , u 2 ) in the triangle determined 
by the equalities. 

0 1 w 2 1 w i , + “ 2 i w 0 / 2 * 
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Number the records j = 1,2»...E. For each of them 
the f oil ow ing operations are performed: 

1) Subtract the average and remove a. linear brand if 
necessary. 

2) Adding zeroes if necessary to obtain a convenient length 
N for the FIT. 

Numbering the observations in the thus modified sample 
from 0 to N-l i x Q x-^ ... x 

3) Fast Fourier transform is performed which gives the 
Fourier coefficients 

N-l 

= l x i exp(-j27ti k/N) 0 ikl N/2 
i=o 

4) Estimate the bispectral density 

k to L w 

(w^, w^-) = , -j-p ) where k,L = 0,l...N/2* 


a) Over a quadratic window 


nf w . w \ A 2 

^ V 2> o 


L L * • 

I l X(o 1 +k 1 )X(w 2 +k 2 )X*(“ 1 + w 2 +k 1 +k 2 ) 

kj^-L k 2 =-l 


or more symmetrically. 


b) over a hoxagona.1 window 


A-i-v ••o = ^ i ?C“i+k 1 )2(“ 2 + k 2> ; 

5M +1 k 1= -l k 2 =-I 

X*( ^ 1 + as 2 + k-^ + kg) 
where the sum is extended over all k^ satisfying 
| k il £1 , I kgl £ L, I l&L+kg I £ 
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5) Finally, average ov-r k pieces 

k 

®("i'“ 2 ) = t 1 YVV 

j=l 

Part 4 is most expensive term in terms of computer time. 

P rocedure A; (Average in the Frequency domain): Evaluate 
operation 5 in the straightforward fashion, e.g. in case (a) 

VVY 4 o- |_j. *Cl +k l) k L X( “2 +k 2 ) 

X*(w 1+ u> 2 + k-^ + k 2 ) 


P rocedure B : (Complex demodulation) 

to 

(1) Put I L (w^) = x(to^ + w k ) for |k| 1 1 

~ o otherwise 


T 

X • L (to k ) = x(to k + <o t ) for IK I £ 21 


k 


0 


otherwise 


This means that we apply a narrow hand pass filter 
and shift frequencies to zero. 

(2) Transform hack- into the time domain. Let n > 2M he 
the least number admitting an efficient FFT. Interpret 

index k modulo n and compute the finite Fourier transforms. 
<o T 03 T 

X ^(s) = I X exp(2^isk/N) 

K to « 

and similarly for X giving x (s). 

3) Average in the time domain: In case (a) compute 
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/\ 


B. (to. , ta T ) A 2 = - 
3 K k L ' o n 


n-1 k w (u +u ) 

l x (s) X (s) X (s) 


S= o 


In case (b) compute 


A 


4M 2 -j n-1 a) oj fto 

j x L (S )'s L (S ) x l (E ) 


3M‘ i +l 


s=o 


It is easy to see that this gives exactly same results 
as the straight forward evaluation of operation 4(a) or ( b ) 
respectively and this is more economical since one does not have 
to compute 'K and 1. 


The third window which is to be discussed is a good 
window. A good window is one which satisfies the following- 
re quire ments ; 

i) The variance of the estimate is as small as possible, 

ii) The width of the window must be as small as possible 

to reduce the computing costs 

iii ) The window must be concentrated near the origin to 
reduce the bias of the estimate. 


The window to be discussed is optimum in the sense 
that the bias of the estimate obtained by using this window 
is the minimum possible [5). The variance of the estimate 
and computing costs are approximately the same in all the 
three cases. This optimum window in the time domain can be 
expressed as follows i 

“ 0 (v V = 'Mu) W Avl 

| | 

where d 0 (x) = ~ | Sin(iu?/M) | + (1 jjp)cos(itT;/M) 
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and the corresponding expression in frequency domain' is as 
follows : 


W n (o 


o^ w l : 


w ) = 
2 ' 


1 

2it 


CO 


V“r 


•y) 


D 0 ( w 2 +y) 


D 0 (y) dy 


whe re D ( 



1+ Cos M 

(M' 2 o 2 — -jt 2 )2 


where M is the maximum time lag of the window. 


This window is used for estimating the hispectral 
density in this work. 


Thus we have discussed the principle behind bispo ctral 
holography and estimation of bispe ctral density in this 
chapter. 
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CHAPTER 4 

COMPUTER SIMULATION OF THE! IMAGING MODEL 

In the previous chapter the methodology for recording 
the hologram using hispectral analysis, and the computational 
aspects involved in the estimation of the bispectra were 
discussed. In this chapter the simulations that have been 
carried out for imaging the objects using bispectral analysis 
are given. 

4 *1 Image Reconstruction of a Point Source 

The model used for computer simulation is shown 
in Fig. 4*1* Here the object considered is a point source. 

The signal emitted by the object is 

5 ( t ) = Cos( ut + <p) +a 2 Cos ( 2 “t + cp) + Cos (4 “ t + 6(t) ) 

(4.1) 

where 6(t) is randomly fluctuated in time with uniform 
distribution over the range (-71 to n). The signal chosen 
above is qua si-periodic, zero-mean and non-Gaussian. The 
noise that is added to the signal at the receivers is assumed 
to be Gaussian. The computer generated pseudo white 
Gaussian signal is taken as noise . The results of simulation 
for the power spectrum and bispectrum of the signal and noise 
are as follows. The power spectrum and bispectrum of both 
the signal and noise are computed. They are shown in 
Fig. 4.2, Fig. 4. 3 and Fig. 4. 4. The Figs. 4-3 and 4.4 are 





o 

o 

o ^ 

o 

c 

o ^ 

O => 
o D~ 


o 

. Q 
jf" UJ 

r~“ , 

(/> < 
z: cr - 

1 1 j u J 
D 

UJ 

• — ^ O i « s 

< 

f 'y- liO 

“•* cc ~~ 

K '“ UJ O 

C ) « * 5 * 

UJ “i 
n r — 
u - Li . z» 

cr Run 
uj uo 


o >* 
tn u 
c 
& 

O D 

o cr 
O <L 
<N w 

Uo 


I 

z 

>-o 

t z 

zy 

UJQ 

DQ 

„jCC„j 

<R< 

U-. «iC 

(Jv't 

UJ<^ 

LjL* **w i ' ' •***»' 

inO~: 

<£ 

rrUJ.- a 

•■*" «-» V I 

Uir’-i/) 

> *r 


OO 



o> 

'-■t 




'I 




« #• * •« 



-4 f! 4 

» o 1 o 


H 

U) 


U) 


CO 

10 

CD 

o 

II 


_ 0> 

D ,JZ 

cr^*^ 

c - £ 
o Q x 
o >o 

W. ID 

CO * 

5 § N 

m on-D x 

o 

>f 1: dO 
O o WU1 

CJ 

xn 'wv-5? c 
rh CD Q. hr 

f *'£ g 

X -X CD CL 


D 

O 

■of-™* 

c 

o 

tJ 

d) 

13 


o 


r-'i 

6 

* iw 

>•' 


j 

CO 


if) 

Cm 


(. 4 # 

CM 


»f> 


CM 


•f> 


t-t t I I I f M i i I I i t i t t * I * t 4 t * » * t 


iil 


"f** 


***• 

<4* 



(\t 

i r . 


0 

t 

<y 

.) 

0 

9 ) 


OJ 


i 



Avi$ua)ui pasipuj jo[\] 




52 


4 • 2 Image Reconstruction of Two Adjacent Point Sources 

The model used for simulation is similar to the 
one shown in Fig. 4-1? except for the fact that? two 
point sources are considered in the object plane instead 
of one. These two point sources are separated toy a distance 
more than the resolution limit 2 xD/x g where D is the 
distance of separation between the object line and scanning 
line. The same procedure as adapted in the case of a single 
point source is followed here. The reconstructed image is 
shown in Fig. 4.5 (b). In this figure reconstructed images 
for two cases are shown. In one case the two point sources 
are assumed to emit identical signals . In the other case 
it is assumed that the signals emitted by the two point 
sources are assumed to have different amplitude distributions. 
It is observed that the spacing between the two point sources 
in the reconstructed image is same as the spacing between 
the two point objects. 

4 • 3 Image Reconstruction of an Uniformly Illuminated - 
Patch 1 

The hologram considered here is a two dimensional 
one. Here also two detectors are employed, one fixed at the 
origin and the other scanning the hologram plane. The set tip 
for this is shown in Fig. 4-6. The signal x(t) emitted by 
the object is expressed as in Eq. (4.1 ) . The contribution 
of the object signal to the signal detected at a point r 
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on the hologram plane is obtained by integrating x( t ) over 
the whole object. As it takes more time to compute the 
integral, to find out the signal contribution the object is 
considered as a matrix (5 x 5) of point sources. Because of 
this fact the reconstructed image has slight intensity 
variations. The reconstructed image . is shown in Pig. 4.7* 

It is obtained by contour mapping. All the contours shown 
in the figure correspond to same height. Since the variations 
in the intensities are very small, the contour of chosen 
intensity appear as shown in the figure. It is not giving 
a clear indication of the object. Because of that Hie Yfr^YUliucted 
vdluei in Tafc>. 4.1 which clearly gives the 

intensity variations. 

In all the above cases, the medium considered causes 
only linear phase shift and attenuation. 

The signal processing steps involved in the data 
acquisition and computer simulation of bispectral imaging 
technique are illustrated in Pig. 4*9* 
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CHAPTER 5 

CONCLUSIONS AND SUGGBSTI ONS 


A study of the paper titled 'Bispectral Holography 1 
hy Takuso Sato and Kimio Sasaki initiated this work. After 
a brief discussion of the principles of holography and 
differences between the optical and acoustical holography > 
the technique of 'bispectral holography' as proposed by Sato 
and Sasaki is examined. A survey of the work carried out in 
this area is done. The fact that' the bispectral ana.lysis can 
eliminate the effects of noise and channel on the signal 
completely' is re-established using computer simulation with 
the help of three examples. In this work the channel is 
assumed to be 'deterministic'. Because of the immediate 
nonavailability of the information on realistic models of 
randomly turbulent media, no example pertaining to this case 
could be included, though the theory applies to this media 
also. This technique naturally requires more computation time 
than the conventional technique as one has to compute both the 
cross- and aut o-bispe ctral densities of the scanned and 
reference signals and their ratio whereas in the conventional 
method one has to compute only the cross-power spectral 
density between the scanned and reference signals, for 
get ting, the hologram signal, but the advantages gained 
clearly justify the increased computation. 
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5 . 1 Suggestions for Future Work 

Apart from its current applications to acoustic 
imaging, passive sonar [63 and machine system diagonosis [73? 
it is felt that the technique can he applied to problems of 
channel characterization and communications. This is 
elaborated in the following discussion. 

a ) Channe 1 Ohara cte rizati on : 

The procedure for channel characterization using the 
present method is as follows. A zero mean, non- Gaussian signal 
is transmitted through the channel to be characterized. 

Assume the corrupting noise to be additive Gaussian and 
statistically independent of the signal. The bispectrum of 
the received signal will be the bispectrum of the channel 
corrupted transmitted signal as the bispectrum of the noise 
becomes zero. The mathematical procedure for finding out the 
transfer function of the channel from the bispectrums of the 
transmitted and received signals is given below. 

Let h(t»'c) be the impulse response of the linear 
time varying channel to be characterized. A known signal 
x(t) is transmitted through the channel. It is assumed to 
be a quasi-peri odic, non- Gaussian zero mean process. Its 
third order correlation function is assumed to be ergodic. 

The bispectral density of x(t) is obtained by observing the 
signal x ( t ) over a period T (T is sufficiently large to 
satisfy the ergodieity property and is given by 
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B X (“V “l> = T x (“k> x ("l> + “i) (5.1) 

The signal received at the other end of the channel 
is given by 

00 

Y(t) = /x(x) h(t,x)dx + n(t) 

— 00 

= m(t) + n(t) 


Assuming that Y(t) is zero mean and its third order 
correlation function is ergodic, R (u.-t ) is expressed as 

y 2 

i T 


T 


I y(t)y(t+T 1 )y(t+T 0 ) (assuming T is sufficiently large 


to satisfy the ergodicity property) 


T 


/m(t) m(t+x-^) rr^t+x^) ('This is true under the assumption 

that m(t) and n(t) are statis- 


tically independent). 

••• B y (» k ,o I| ) = | M(» k ) M^) M*(» k +U. L ) 

= i Z(“ k ) H(» k ,t) !(%) H(» I( t) X> k + L ) 

H(“ k .“i,t) • (5.2) 

(Assuming that the channel is not varying over the interval 
of observation) where t is any time instant in the interval T« 
The ratio of Bets. ( 5.1) and (5*2) gives 

B h (<V ( V t)=H( ‘V X ) H(o. l ,x) H*(o fe +u) L> x) 

Evaluate B H (o k ,u L ,x) for k»L = 1,2, ...n. 

Then arranging these values in matrix form 



... t)H(w nr T)H(i*i n +w- L ) t) 

H(» 2 ,t)H(» 2 ,t)H*(» g+“g. t) 


V T)H <“n' T > H ^“n + V' t >_ 

> x > 

B^VVO _ 

To solve for £[( 1 ^, 1 ;), H(u 2 ,t) H(o 2n j.T) there are only 

n ( ii+l ) 

~ 2 — ^ equations since the above matrix is a symmetric 
matrix. 

Since the above system is an over determined system 
of non-linear equations with multiplicative terms to find 
a solution one will have to firt linearize them by taking 
the logarithm and then resorting to any one of the methods 
of solution for over determined system, e.g. pseudo inverse 
[16) » least squares or residual projection method [ 173 * 

b) Communication through a Bandomly Time Varying Media: 

The image of the object re constructed using bispe ctral 
analysis has the effect of neither the channel nor the additive 


H(ai n ,T)H(w 1 ,T)H*(a) n -Ho 1 ,T) H( 

* x ) W 1 * U 2 ' T ) 


B tt ( w ,0 , t ) 
L H v n* n J ' 
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Gaussian noise. The bi spectral analysis may be applied in 
the communications contest where the aim is to obtain the 
signal free from additive noise and the effects of the 
channel. To get the ima.ge of the object the variation 
of the amplitude and phase of the object wave relative to 
those of the reference wave in space has to be determined 
whereas to get the signal transmitted, the variation of the 
amplitude and phase of the object signal with time has to 
be determined. Por obtaining the signal, an analysis 
similar to the one used for getting image of the object 
ca.n possibly be carried out by interchanging time and space 
domains . 



APPENDIX 


This contains some of the important listings 
of the computer simulation. Comment cards introduced 
in ach routine explains the function of that routine. 
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